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Orientation-dependent binding energy of graphene on palladium 
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Using density functional theory calculations, we show that the binding strength of a graphene monolayer on 
Pd(lll) can vary between physisorption and chemisorption depending on its orientation. By studying the 
intcrfacial charge transfer, we have identified a specific four-atom carbon cluster that is responsible for the 
local bonding of graphene to Pd(lll). The areal density of such clusters varies with the in-plane orientation 
of graphene, causing the binding energy to change accordingly. Similar investigations can also apply to 
other metal substrates, and suggests that physical, chemical, and mechanical properties of graphene may be 
controlled by changing its orientation. 



Since its isolation in 2004, 1 graphene has attracted 
a great deal of attention because of its unique phys- 
ical, chemical, and electronic properties.- 2 ' - " A remain- 
ing technological hurdle that prevents the introduction 
of graphene into specific applications is a basic, pre- 
dictable understanding of the interaction of graphene 
with supporting substrates £ The sensitivity of graphene 
to the substrate material and surface orientation is well 
established, but even relative azimuthal orientation of 
graphene on these substrates can affect its properties^— 
These recent works address manifestations of the interac- 
tions between graphene and its host substrate, showing 
that these interactions actually vary with the local envi- 
ronment, that is, the chemical or electronic environment 
in the immediate vicinity of certain carbon atom(s). 11 
Careful bookkeeping of the local structures (and the asso- 
ciated interactions) present within a single spatial period 
of a (moire) superstructure can be compared to varia- 
tions in the global properties, such as binding energy, 
to deduce the impact of the orientation of the graphene 
layer. 

In this letter, we study monolayer graphene on Pd(lll) 
for a range of in-plane orientations, and show that the 
atomic-level carbon-palladium registry affects the bind- 
ing to the substrate through changes to molecular or- 
bitals at the graphene-metal interface. Furthermore, we 
uncover a link between the interfacial binding strength 
and the areal density of certain four-atom clusters cen- 
tered atop Pd atoms in the first substrate layer. This 
link provides fundamental insight into the physical ori- 
gin of the orientation dependence of the binding energy, 
and can be exploited for understanding the interaction of 
graphene with other substrates as well. 

In order to study the orientation dependence of the 
binding energy, we start by constructing perfectly peri- 
odic moire patterns of graphene on Pd(lll). All moire 
patterns are incommensurate when using the equilibrium 
lattice constants of graphene and Pd, so we have en- 



sured commensurability at each orientation (defined as 
in Fig. [T]) by applying small strains to the graphene lat- 
tice. We have performed density functional theory (DFT) 
calculations of graphene domains having orientations 
between and 30° with respect to the substrate; other 
angular domains reduce to < 8 < 30° through rotational 
symmetry. For convenience, we have chose four different 
orientations, — 5.7, 10.9, 19.1, and 30.0°, the first three 
of which have been experimentally identified in Ref. [T2I ; 
the 30° orientation helps compare our results with those 
of Giovannetti et al. e To understand the impact of atomic 
registry on the graphene-Pd(lll) binding, each carbon 
atom is designated as a top (T), bridge (B), gap (G), 
or hollow (H) site, as shown in Figure [TJ based on its 
position relative to the substrate. 
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FIG. 1. (Color online) Possible occupancy sites for carbon 
atoms on Pd(lll): T (top, C atom atop a Pd atom), B 
(bridge, C atom directly above the midpoint of a surface Pd- 
Pd bond), H (hollow, C atom above the hollow site), and G 
(gap, C atom at the center of the smallest TBB triangles). 



Total energy calculations were performed using the 
Vienna Ab-initio Simulation Packag e 13 ' 14 with ultrasoft 
pseudopotentials in the local density approximation, 15 
with a 286 eV plane- wave energy cutoff, and a 12 A vac- 
uum thickness. The Brillouin zone was sampled using a 
Monkhorst-Pack grid-7 x 7 x 1 (5.7°), 9x9x1 (10.9°), 
17 x 17 x 1 (19.1°), and 35 x 35 x 1 (30.0°)-with more 



2 



than 100 points/A" 1 along each reciprocal lattice vec- 
tor. The substrate was modeled with four Pd(lll) layers 
(lowest two kept fixed during relaxation) , with the excep- 
tion of the 5.7° system, which was modeled with three Pd 
layers due to computational limitations; we have tested 
that binding energy trends obtained with 3 or 4 substrate 
layers are indeed consistent with one another. We have 
computed the binding energy per unit area E^ via 



VPd 
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GrPd 
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where i?GrVPd is the total energy of a system in which 
the graphene sheet is far from the Pd surface, and i?GrPd 
is the total energy of a system where the graphene layer 
was relaxed on the substrate. We express the area A in 
carbon atoms rather than A 2 (1 C = v / 3a 2 /4, where a is 
the lattice constant of graphene). 

Binding energy is a fundamental property of epitax- 
ial graphene systems, affecting not only its structural 
and mechanical stability, but the electrochemical prop- 
erties as well. While calculations of binding energy have 
become commonplace, our goal here is to correlate it 
with simple geometric features that are present in the 
moire superstructures characteristic to various orienta- 
tions. Seeking such correlation involves, in order: (a) 
identifying the possible locations (coincident sites) that 
a carbon atom in graphene can occupy with respect to 
the substrate, (b) tracking the populations of these co- 
incident sites (i.e., their areal densities), or of groups 
of them, as functions of the orientation of the graphene 
sheet, and (c) finding which of these tracked popula- 
tions depend on the orientation angle in a manner sim- 
ilar to the binding energy. Following this plan, we start 
by selecting four distinct types of coincident sites to de- 
scribe the positions of a carbon atom relative to the Pd 
surface-top (T), bridge (B), hollow (H), and gap (G), 
as described in Figure [1] To avoid counting ambiguities 
when tracking site populations, we have chosen to iden- 
tify the regions for which coincident sites are populated 
by site-centered circles that touch without overlapping, 
as shown in Fig. [Q. This convention ensures that 90.7% 
of the substrate surface is counted; atoms within the re- 
maining 9.3% of the surface (colored dark blue in Fig. [1]) 
are disregarded. 

Tracking changes in coincident site population requires 
a fixed reference, and we have selected this reference to 
be the expected fraction of carbon atoms that would lie 
at a certain coincident site under the assumption of ran- 
dom distribution of carbon atoms. Placed at random, 
1/12 of carbon atoms lie on top sites and 1/4, 1/6, and 
1/2 at bridge, hole, and gap sites, respectively. Figure [5] 
shows the variation in the population of each coincident 
site as a function of angle. Although a given population 
of single-type coincident site may vary by as much as 45% 
with respect to its reference, no clear trends emerge for 
the variation of any coincident site population with the 
orientation 8 (Fig. [2J. On the other hand, as we show 
below, the binding energy of graphene monolayer on Pd 
varies monotonically with the orientation angle: there- 
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FIG. 2. (Color online) The areal density of the individual 
types of sites occupied by the C atoms of the graphene sheet 
as a function of orientation. 



fore, no correlation exists with between binding energy 
and populations of any of the single coincident sites, T, 
B, H, or G. 
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FIG. 3. (Color online) Graphene on Pd(lll) oriented at 
9 = 5.7°, 10.9°, 19.1°, and 30.0°, with the carbon atoms 
colored according to the site they occupy (defined in Fig. [T|). 
The moire surface unit cells are outlined in white. Despite the 
nearly random distribution of top, bridge, hole, and gap sites, 
we note the periodic repetition of distinctive clusters of sites 
such as TB3 [three bridge site carbons (yellow) surrounding 
a top site (red)], TH3 (three green, central red), BG3 (three 
cyan, central yellow), HT3 (three red, central green), HH3 
(three green, central green), HG3 (three cyan, central green), 
GB3 (three yellow, central cyan), and GG3 (three cyan, cen- 
tral cyan). 



A simple inspection of the moire patterns correspond- 
ing to the orientations 8 studied here (Figure [3]) suggests 
that coincident sites are occupied in a coordinated man- 
ner, which is expected since the carbon atoms that oc- 
cupy them are part of the same graphene lattice (refer to 
examples given in the caption to Fig. . Therefore, we 
proceed to consider nearest neighbors of a given (central) 
site, along with that site, in order to create an inventory 
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FIG. 4. (Color online) (a) Density of different types of 4-site 
clusters for the orientations of 9 = 5.7°, 10.9°, 19.1°, and 
30.0°. (b) The binding energy Et of graphene on Pd(lll) in- 
creases as a function of orientation 9. The only cluster whose 
areal density increases monotonously with 9 is TB3, which 
we show to be responsible for the bonding of graphene to the 
substrate. 
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ure [5] shows the p z states of the carbon atoms in a TB3 
cluster and the d states of the Pd atom beneath this clus- 
ter. The presence of common peaks for these projected 
densities of states (marked by vertical gray lines in Fig. [5]) 
below the Fermi level and close to it indicates the forma- 
tion of hybridized, bonding orbitals between the p z states 
of the carbon atoms of the cluster and the d-states of the 
Pd atom underneath. 
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FIG. 5. (Color online) Site and angular momentum projected 
density of states (PDOS) for the top and the bridge site car- 
bons in a TB3 cluster and for the corresponding Pd atom at 
the 30.0° orientation. The presence of common peaks (verti- 
cal gray lines) for the C-p z , Pd-d xz , Pd-d^z, and Pd-d z 2_ r 2 
projected density of states is consistent with the formation of 
hybridized orbitals at the interface. 



of clusters that can occur on the substrate. We have 
identified these clusters by a two-letter abbreviation; for 
example, HG3 is a hollow-site (H) carbon surrounded by 
three gap-site (G3) nearest neighbors. Eight such clus- 
ters are represented in the moire patterns shown in Fig.|3l 
and we have determined the areal density as a function 
of angle 6 for all of them [Fig.HJa)]. As seen in Fig.rjjb), 
only one of these four-atom clusters, TB3, has a density 
that varies with the orientation angle in the same way 
in which the binding energy does. At small orientations, 
the graphene sheet physisorbs to the Pd(lll) substrate 
with a binding energy of only 41 meV/C. As the sheet is 
rotated and the number of TB 3 clusters increases, so does 
the binding energy, increasing to 73 meV/C at 6 — 30° 
where the density of TB3 clusters reaches its maximum 
of 9.65 x 1CT 2 A" 2 . 

To substantiate the link between areal density of TB3 
clusters and binding energy beyond the comparison of- 
fered by Fig. Q|b), we have analyzed the site-projected 
density of states (PDOS) corresponding to the atoms of 
the TB3 cluster, as well as the electronic transfer oc- 
curring upon the creation of the graphene-Pd interface. 
Given the need for very fine Brillouin-zone sampling, the 
PDOS calculations were performed on the relaxed struc- 
tures with the order- N code SIESTA 16 using a double-^ 
basis set and a 70 x 70 x 1 Monkhorst-Pack grid. Fig- 



Electronic transfer calculations also reveal a clear sig- 
nature of the bonding between graphene and Pd for a 
certain range of the orientations. At small angles (6 < 
10.9°), there is no significant charge transfer, consistent 
with the low binding energy values shown in Fig. IDJb). 
However, for 9 = 19.1° and 30° the electron transfer be- 
comes significant, as shown in Figs. [6^a) and (c): this 
transfer amounts to the formation of chemical bonds 
(occupied bonding orbitals) between graphene and sub- 
strate, which are the physical origin of the binding energy 
increase computed at large angles [Fig. BJb)]. In Fig- 
ure EJa), we can identify three bonding orbitals by their 
shapes: the axially-symmetric TB3 bond, the ovoid bond 
that lies below the cr-bond between two carbon atoms, 
and the oblong bond where adjacent carbon sites are not 
directly atop a first-layer Pd atom. The "strengths" of 
the bonds formed, as estimated by the increased elec- 
tronic charge, Q, in a certain volume,— are also dif- 
ferent with the highest corresponding to the TB3 case: 
Qtb 3 = 0.055e > Q 

ovoid — 

0.042e > Qobiong = 0.007e. 
This reinforces, albeit qualitatively, our earlier observa- 
tion that the combined effect of the top and bridge-site 
neighboring carbons is responsible for bonding. Unlike 
the case of the 19.1° system where three types of bonds 
are found, in the case of 30° only TB3 is present [Fig- 
ure E^b)]. The TB 3 bonds are ~24% weaker for 6 = 30° 
than their 19.1° counterparts, but are present on 2/3 of 
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the surface Pd atoms which accounts for the the higher 
binding energy. 




• • • 

ovoid oblong TB3 

FIG. 6. (Color online) Electron transfer density for (a) 19.1° 
and (b) 30.0° orientations, showing regions of charge accumu- 
lation (red-orange spectrum) and depletion (blue spectrum). 
Three types of bonds can be identified at these orientations, 
marked in (a) as ovoid, oblong, and TB3. Only the TB3 
type exists at 9 = 30° (b), and no bonds are formed at orien- 
tations of 5.7° or 10.9°. Panel (c) shows side-views of the 
charge transfer regions corresponding to different types of 
bonds identified in (a). 



When bound strongly to a substrate, graphene tends 
to adopt one predominant orientation, as seen on 
Ru(0001)^ and Ni(lll). 7 On substrates with weaker 
binding, including Pd(l 1 1) P- Pt (1 1 1 ) and Ir ( 1 1 1 ) ,21 
multiple azimuthal orientations are observed, and are 
readily identified by the presence of moire structures 
with different spatial periodicities L. Recently, Merino et 
al££ have catalogued structures of graphene on Pt(lll) 
into minimum-strain "phases" , characterized by their 
orientation-dependent spatial periodicities. It is worth 
noting that when applied to the graphene on Pd(lll) 
system, the Merino et al. model describes well the struc- 
tures we examined here: 5.7°, L = 17.12 A (same as in 
Ref.lH C phase); 10.9°, L = 11.20 A (10.75 A in Ref.[H, 
k phase); 19.1°, L = 7.33 A (7.25A in Ref. 22, (i phase); 
and 30°, L = 4.89 A (5.00A in Ref. [H, a phase). The 



most stable phases (orientations) are those with stronger 
binding, i.e., a and /3, as shown before in Fig. SJb). 

In summary, we have shown that the binding energy of 
graphene on Pd(lll) depends on its orientation, and that 
this dependence arises not from the population of any 
single coincident site, but from coincident site four-atom 
clusters (one top-site C atom surrounded by three bridge- 
site carbons). The 78% increase in binding energy, from 
41 meV/C to 73 meV/C, emphasizes the sensitivity of 
this property to the carbon-palladium coincidence. The 
approach presented here can be applied to other metal 
substrates as well, in order to see if equally simple clus- 
ters (and which ones) are responsible for controlling the 
binding strength between physisorption and chemisorp- 
tion. 
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